Stochastic Analysis and Regeneration of Rough Surfaces 
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We investigate Markov property of rough surfaces. Using stodiastic analysis we characterize the com- 
plexity of the surface roughness by means of a Fokker-Planck or Langevin equation. The obtained 
Langevin equation enables us to regenerate surfaces with similar statistical properties compared 
with the observed morphology by atomic force microscopy. 



Studying the growth, formation and morphology of in- 
terfaces has been one of the recent interesting fields of 
study because of its high technical and rich theoretical 
advantages [1]. One of the main problems in this area is 
the scaling behaviour of the moments of height difference 
Ah = h{xi) — h{x2) and the evolution of the probability 
density function (PDF) of A/i, i.e. P{Ah, Ax) in terms of 
the length scale Ax. Recently Friedrich and Peinke have 
been able to obtain a Fokker-Planck equation describing 
the evolution of the probability distribution function in 
terms of the length scale, by analyzing some stochastic 
phenomena, such as turbulent free jet, etc. [2-4]. They 
noticed that the conditional probability density of field 
increments (velocity field, etc.) satisfies the Chapman- 
Kolmogorov equation. Mathematically this is a necessary 
condition for the fluctuating data to be a Markovian pro- 
cess in the length scales [5]. 

In this letter using the method proposed by Friedrich 
and Peinke, we measure the Kramers-Moyal's (KM) co- 
efficients for the fluctuating fields Ah and h{x) of a de- 
posited copper film. It is shown that the first and sec- 
ond KM's coefficients have well-defined values, while the 
third and fourth order coefficients tend to zero. There- 
fore, by addressing the implications dictated by the the- 
orem [5] a Fokker-Planck evolution operator has been 
found. The Fokker-Planck equation for P{Ah, Ax) is 
used to give information on changing the shape of PDF 
as a function of the length scale Ax. By using this 
strategy the information of the observed intermittency 
of the height fluctuation is verified [6]. The first and 
second KM's coefficients for the fluctuations of h{x), en- 
ables us to write a Langevin equation for the evolution 
of height with respect to x. Using this equation we re- 
generate the surface with similar statistical properties, 
compared with the observed morphology by atomic force 
microscopy. The regeneration of a surface is known as 
the inverse method. There are other inverse method ap- 
proaches introduced in the literature [13]. In the previous 



attempts, to regenerate the surface, an evolution equa- 
tion for h{x, t) vs t has been evaluated. Here we do this 
by an evolution equation for h(x) vs x, for a certain time. 

For this purpose, a copper film was deposited on 
a pofished Si(lOO) substrate by the resistive evapora- 
tion method in a high vacuum chamber. The pres- 
sure during evaporation was 10"^ Torr. The thickness 
of the growing films was measured in situ by a quartz 
crystal thickness monitor. We performed all deposi- 
tions at room temperature, with a deposition rate about 
20 — iQnm/min. The substrate temperature was deter- 
mined using a chromel/alumel thermocouple mounted in 
close proximity of samples. The surface topography of 
the films was investigated using Park Scientific Instru- 
ments model Autoprobe CP. The images were collected 
in a constant force mode and digitized into 256 x 256 pix- 
els with scanning frequency of 0.6 Hz. The cantilever of 
0.05 N spring constant with a commercial standard 
pyramidal Si^N^ tips was used. A variety of scans, each 
with size L were recorded at random locations on the Cu 
film surface. 

It is a common procedure to characterize the complex- 
ity of a rough surface by checking the scaling behaviour 
of the moments C<j =< \h{xi) — h{x2)\'^ > in terms of the 
length scale Ax = |xi — X2|. We investigated the scal- 
ing behaviour of the q-th moment Cq and observed that 
all of the moments (up to q = 20) behave as |xi — X2|^'' 
within the scaling region ^ 10 to 150 nm. We have found 
a nonlinear relation between S,q and q. This shows that 
the height fluctuations are intermittent or multi-fractal 
( see [8, 15 ] and references therein). The roughness 
exponent a is related to the exponent ^2 as a = ^2/2 
[1]. For the stationary samples with thickness 440 nm, 
the roughness exponent a was found to be 0.83 ± 0.03. 
From the stochastic point of view one has to remark that 
multi-fractality is based on properties of the roughness on 
distinct length scales. However, checking the scaling be- 
haviour does not explain possible correlation between the 
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FIG. 1. Test of Chapman-Kolmogorov equation for differ- 
ent values A/ii = — 21nm, A/ii = and A/ii = 21nm. The 
bold and open symbols represent directly evaluated PDF and 
the integrated PDF, respectively. The length scales Axi, Ax2 
and Ax3 are 180nm, 320nm and 260nm, respectively. 

roughness measures on different scales. Also it is noted 
that the methods based on multi-fractality are limited to 
the subclass of rough surfaces which show scaling prop- 
erties. The method introduced by Friedrich and Peinke 
is a general method, which explains the complexity of 
the surface roughness, with no scaling feature to be ex- 
plicitly required. Their method yields an estimation of 
an effective stochastic equation in the form of a Fokker- 
Planck equation (also known as Kolmogorov equation). 
The connection between the multifractality and Marko- 
vianity has been discussed in [4]. 

A complete characterization of the statistical proper- 
ties of the height fluctuation requires the evaluation of 
joint PDF's PAr(A/ii, Axi; Aft,jv, Accjv), for any ar- 
bitrarily N. If the process is a Markov process (a pro- 
cess without memory), an important simplification arises. 
For this type process the N-point joint PDF, Pjy, is 
generated by a product of the conditional probabilities 
P(A/i,+i,Aa;,+i|A/i,,Aa;,), fori ^ 1,...,7V-1. As 
a necessary condition for being a Markov process, the 
Chapman-Kolmogorov equation, 

p(A/i2, Ax2|A/ii, Axi) = 

J d{Ahs)p{Ah2,Ax2\Ah3,Ax3)p{Ahs,Ax3\Ahi,Axi) (1 

should hold for any value of Ax^, in the interval 
Ax2 < Axa < Axi [5]. We checked the validity of the 
Chapman-Kolmogorov equation for different Ahi triplets 
by comparing the directly evaluated conditional proba- 
bility distributions p{Ah2, Ax2\^hi, Axi) with the ones 
calculated according to rhs. of eq.(l). In Fig. (1), the 
two direct and integrated PDF's are superimposed for the 
purpose of illustration. The bold and open symbols rep- 
resent directly evaluated PDF and the integrated PDF, 
respectively. Assuming a statistical error of the square 
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FIG. 2. Drift and diffusion coefficients D^{Ah) and 
D'^{Ah) are estimated from the eq. (3). The and 
present the linear and quadratic behavior, respectively. 

root of the number of events of each bin we find that 
both PDF's are statistically identical (see also [14] an 
another interesting and carefully presented example of 
application of the Chapman-Kolmogorov equation). 

It is well-known, the Chapman-Kolmogorov equation 
yields an evolution equation for the change of the dis- 
tribution function p{Ah, Ax) across the scales Ax. The 
Chapman-Kolmogorov equation formulated in differen- 
tial form yields a master equation, which can take the 
form of a Fokker-Planck equation [5] : 

["MX^^'^^^'^'^) + ^D(^H^h,rMAh,r) (2) 

where r := Ax. The drift and diffusion coefficients 
D^^^ (Aft,, r), Z?(^)(Aft,,r) can be estimated directly from 
the data and the moments M^''^^ of the conditional prob- 
ability distributions: 

D'^''\Ah,r) = ^ lim^r^oM^''^ 

Af f*^' = ^ I dAh'{Ah' -Ahfp{Ah',r + Ar\Ah,r). (3) 
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FIG. 3. Probability densities of the height difference 
A/i = h{x + Aa;) - h{x) for the length scales Aa; = 200, 300 
and 500 nm (from top to bottom). The results obtains from 
the data analysis of the AFM image (R) and numerical inte- 
gration of an effective Fokker-Planck equation (E) i.e. eq.(2), 
respectively. The PDF's are shifted in vertical directions for 
convenience of presentation and A/i's are measured in units 
of the standard deviation of Ah at Aa; = 200nm. 



The coefficients D^^\l\h,rys are known as Kramers- 
Moyal coefficients. The drift and diffiision coefficients 
Z)'^-' and Z?*-^-* are displayed in fig. (2). It turns out that 
the drift term Z?^^) is a linear function of A/i, whereas 
the diffusion term D^"^^ is a function quadratic in Aft,. 
For large values of Aft, our estimation becomes poor and 
thus uncertainty increases. From the analysis of the data 
set we obtain the following approximation; 

L)(^^(Aft,Aa:) = -0.0055 Aft 

Li(2)(Aft,Aa;) = [(2.9 x IQ-^Mif + Qmb(^x)^-^^]/ Ax 

where Aft is measured in units of the standard devi- 
ation of Aft at Ax = 200nm. According to Pawula's 
theorem, the Kramers-Moyal expansion stops after the 
second term, provided that the fourth order coefficient 
D^^^ (Aft, Aa;) vanishes [5]. The forth order coefficients 
D^^^ in our analysis was found to be about D*-^-* ~ 
lO^'*!?^^-'. In this approximation we can ignore the coef- 
ficients for n > 3. To perform a quantitative test 
of the result with these coefficients, we solve the Fokker- 
Planck equation for the PDF at scales Ax <C L with a 
given distribution at sample size L [6,7]. Fig. (3) shows a 
comparison between the analysis of AFM image and the 
solutions of the obtained Fokker-Planck equation for the 
copper surface for the length scales Aa; = 200, 300 and 
500 nm. The figure shows that the solutions of our model 
fit the experimentally determined PDF's with good pre- 
cision. In the integral scale our measured PDF is nearly a 
Gaussian distribution. In our approximation the stochas- 
tic process underlying the height fluctuation changes is a 
linear stochastic process with multiplicative noise. 

By the same procedure, we checked the Markovian 



nature of the fluctuations of the height h = ft(a;) — ft, 
and found the following expression for the £1(1' (ft) and 
£1(2) (/i): 



-0.01ft 



0.1 



0.004ft + 5.19 X lO-'^ft^. 



(5) 



The height field is measured in units of the standard 
deviation of ft. Analogous to equation (2), we can write 
a Fokker-Planck equation for the PDF of ft by replac- 
ing r and Aft with x and ft, respectively. We note that 
this Fokker-Planck equation is equivalent to the following 
Langevin equation (using the Ito interpretation) [5] : 



^ft(a;) = Z?(i) (ft) + y^Z?(2)(ft) f{x) 



(6) 



Here, f{x) is a random force, zero mean with gaussian 
statistics, (5-correlated in a;, i.e. < f{x)f{x') >= S{x — 
x'). Furthermore, with this last expression, it becomes 
clear that we are able to separate the deterministic and 
the noisy components of the surface height fluctuations 
in terms of the coefficients D'^^^ and D^^). Equation (6) 
enables us to regenerate rough surfaces which are similar 
to the original one (in the statistical sense ). In fig. (4) 
the AFM and regenerated images are demonstrated. The 
regenerated surface is very similar in statistical sense to 
the original one. All regenerated patterns are statistically 
similar. To ensure this fact, for instance, in Fig. (5) we 
have plotted the second moment of the structure function 
6*2 for the AFM and a regenerated surfaces and their 
roughness exponents were found 0.83 ± 0.03 and 0.83 ± 
0.01, respectively. 
(4) There are a few comments on the regeneration of rough 
surface that we would like to notify. When we are dis- 
cussing about a Markov process, one should note that, 
this is true within an approximation. For instance, in 
random motion of a particle inside a fluid, it is known 
that the collision of the particle with the fluid molecules 
is not instantaneous, and takes a certain duration. Dur- 
ing the time that a collision is taking place, the change of 
velocity is not Markov, because the velocities in the col- 
lision time scale have memory. Consequently, in the time 
series for the velocities of the particle, if the time inter- 
vals are less than the collision time scale, the process can 
not be regarded as a Markov process. The minimum time 
interval that the particle motion can be considered as a 
Markov process is known as Markov time scale and the 
motion is known as a Brownian motion. In the stochas- 
tic analysis of the rough surface, we are dealing with the 
Markov property of height fluctuations in spatial dimen- 
sions, therefore, instead of a Markov time scale, here we 
will have a Markov length scale Imarkov Our analysis 
shows that Imarkov — 160 nui, equivalent to 8 pixels in 
our AFM image [9]. The surface is regenerated by iter- 
ating the eq.(6), which gives us a series of data without 
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FIG. 4. AFM and regenerated surface images (from up 
to bottom), which we have regenerated the rough sur- 
face using the Langevin equation for dynamics of h(x). 
As drift term D^^^ (h) = — O.Ol/i and as diffusion term 
D'^'(/i) = 0.088 - 0.004ft + 5.19 x 10"®/i^ 



memory. To compare the regenerated surface with the 
original one. we have to take the spatial interval in the 
numerical discritization of eq.(6); to be equal to 1 pixel. 
However, here the Markov length is equal to 8 pixels. 
Therefore, we should relate the height field within the 
Markov length. There are a number of methods to corre- 
late the generated data in this interval [9] . We do this by 
means of scanning the surface with a tip, where its size is 
about the Markov length [12]. The tip that we have used 
has the form z = ax^ + by^. In this case the parameters 
a and b are 0.035. 

In summary, we have shown that the probability den- 
sity of height increments satisfy a Fokker-Planck equation 
which encodes the Markovian property of these fluctua- 
tions in a necessary way. We are able to give the ex- 
pression of Kramers-Moyal's coefficients for the stochas- 
tic processes Ah and h by using the polynomial ansatz 
[10,11]. Also we could find the form of path probabil- 
ity functional of the height increments in spatial scale, 
which naturally encodes the scale dependence of proba- 
bility density. This gives a clear picture about the inter- 
mittent nature in height fluctuations. The methods en- 
ables us to regenerate many realizations of the rough sur- 
face with similar statistical properties in favored scales. 
As an application, large surface generation would be pos- 
sible by sampling the real surface with high resolution (in 
same the resolution as nanoscope imaging, e.g. AFM im- 
ages). This would be applicable in computer simulation 
of the surface and interface processes, for example, the 
diffusion of materials between rough surfaces, the effect 
of roughness on the friction, and so on. 

We thank F. Azami, M. Sahimi, and N. Taghavinia for 
their useful discussions. 
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FIG. 5. Log-log plot of the second moment of 
height-difference vs I, for real sample and regenerated sam- 
ple. The roughness exponents for real and regenerated are 
0.83 ± 0.03 and 0.83 ± 0.01, respectively. 
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